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Abstract 

In quasi-single field inflation models, massive isocurvature modes, that are coupled to 
the inflaton and have mass of order the Hubble parameter, can have nontrivial impacts on 
density perturbations, especially non-Gaussianities. We study a simple example of quasi- 
single fleld inflation in terms of turning inflaton trajectory. Large bispectra with a one- 
parameter family of novel shapes arise, lying between the well-known local and equilateral 
shape. The trispectra can also be very large and its magnitude tjvL can be much larger than 
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1 Introduction and summary 



Inflation [Ij has become the leading paradigm of the early universe. However, the detailed 
dynamics of inflation is still a mystery. A major theme in cosmology is to build inflationary 
models and compare their predictions with experimental data. The theoretical predictions 
of general single field inflation models have been well understood. Nonetheless, there exists 
another important possibility that the inflationary dynamics can involve multiple flelds. This 
leads to a variety of new models and interesting phenomenologies. 

When multiple flelds are involved, the fleld space can be decomposed into the inflationary 
direction and isocurvature directions. The quanta in these fleld directions are called inflaton 
and isocurvatons, respectively. In this paper, we investigate a class of models where there 
is one flat slow-roll direction, and all the other isocurvature directions have mass at least 
of order the Hubble parameter H. We call this class of models quasi-single field inflation. 
If the inflaton decouples from the isocurvatons or the isocurvaton mass are all much larger 
than 0{H), quasi-single fleld inflation makes the same prediction as the single fleld inflation. 
However, once large couplings exist and the mass are of order H, we will show that these 
massive isocurvatons can have important effects on density perturbations. 

In this paper, we shall study a simple model of quasi-single fleld inflation [2]. In this 
model, the coupling between the inflaton and the massive isocurvaton is introduced by a 
turning trajectory. The tangential direction of this turning trajectory is the usual slow-roll 
direction, while the orthogonal direction is lifted by a mass of order H. See Fig. [1] 




Figure 1: This flgure illustrates a model of quasi-single fleld inflation in terms of turning 
trajectory. The 6 direction is the inflationary direction, with a slow-roll potential. The a 
direction denotes the isocurvature direction, which typically has mass of order H. 

The motivations for investigating quasi-single fleld inflation follows. 
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• UV completion and fine-tuning in inflation models. To satisfy the conditions for in- 
flation, fine-tunings or symmetries should generally be evoked. At least this is found to be 
the case for models that have reasonable UV completion in string theory and supergrav- 
ity [3111]. For slow-roll inflation, this means that, in the inflationary background, the light 
fields will typically acquire mass of order the Hubble parameter H, which is too heavy to be 
the inflaton candidates. On the other hand, in a UV completed theory, multiple light fields 
arise naturally. Taking these facts into consideration, a natural picture of inflation emerges: 
There is one inflation direction with mass m ^ if, and some other directions in the fleld 
space with m ~ if. In contrast to the slow-roll potential, large higher order terms in the 
potential such as V" ~ H and V"" ~ 1 can arise naturally in these non-flat directions. To 
have more than one flat direction needs extra flne-tuning. The above picture for inflation 
suggests the quasi-single fleld inflationary models. 

• Inflationary phenomenology. When the inflaton trajectory turns in the fleld space, 
isocurvature perturbation is converted to curvature perturbation. Precisely understanding 
the effect of such a conversion on density perturbations is in general an important but difficult 
question. Quasi-single field infiation provides explicit examples where such a conversion can 
be calculated from first principles and give non-trivial predictions, such as new shapes of 
large non-Gaussianities and running of the density perturbations. 

• Filling the gap between single and multiple field inflation. Quasi-single field infiation 
fills the gap between single field infiation and multi-field infiation, with new observational 
consequences. The relation between single field infiation, quasi-single field infiation and 
multiple field infiation is illustrated in Fig. |2j 

• Non-Gaussianity. Non-Gaussianity is potentially one of the most promising probes of 
infiation. Generally speaking, observing shapes and running of non-Gaussianities can tell 
us about the details of infiatons and isocurvatons during infiation. In the most well-studied 
models, scale-invariant non-Gaussianities usually take either local shape or equilateral shape. 
In quasi-single field infiation, large non-Gaussianities with a family of new shapes lying be- 
tween the local and equilateral shape can arise naturally, and they are potentially observable 
in the near future. This provides new opportunities for interactions between model building 
and data analyses. 

• Methodology. We compute the correlation functions of the curvature perturbation using 
the in- in formalism P,[6]. The methodology used here is inspiring in several aspects. Firstly, 
we use the transfer vertex to account for the transformation from the isocurvature mode to 
curvature mode, and use Feynman diagrams to compute the correlation functions perturba- 
tively. As a result, the leading order contribution comes from the fourth (for bispectrum) or 
higher (for trispectrum) order perturbation theory. Such methods can have broader applica- 
tions in, for example, multi-field infiation models. Secondly, we will use two representations 
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Figure 2: The relationship between single field inflation, quasi-single field inflation and multi- 
field infiation. Here we take inflation evolving two fields as an example. Quasi-single field 
infiation fills the gap between single field infiation and multiple field infiation. 

of the in-in formalism, i.e. the factorized form and commutator form. We find that each has 
its computational advantages and disadvantages and they are complementary to each other. 
Based on this, we develop a "mixed form" of the in-in formalism, which has the advantages 
of both forms. 

This paper is organized as follows. 

In Sec. [21 we study in detail an explicit model of quasi-single field infiation with one 
isocurvature direction, proposed in p]. We investigate the zero-mode solution, the pertur- 
bation theory, solve the mode functions and discuss the transfer vertex (as illustrated in 
Fig. M^a)) between the infiaton and isocurvaton. 

In Sec. [3l we study the power spectrum of this model. The power spectrum receives 
correction from two transfer vertices, as illustrated in Fig. El^b). In the slow-turn case, we 
shall show that the correction is of order 



In Sec. m we calculate the bispectra of quasi-single field infiation. Since the slow-roU 
condition is not required in the isocurvature direction, there is a sharp contrast between 
V" for the massive isocurvaton potential and V^" for the slow-roll potential. The slow-roll 
condition requires Vs/" ~ 0{e^)Pl^'^ H , (where e collectively denotes slow-roll parameters). 




(1.1) 



where is the power spectrum, 6/H is the turning angular velocity in Hubble time. 
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Figure 3: Feynman diagrams for the transfer vertex (a), corrections to the power spectrum 
from isocurvature modes (b), and the leading bispectrum (c). 



which contributes an 0{e^) term to f^L [Z]- However V'" is not subject to such conditions. 
For example, V" can naturally become of order H. To estimate the size of the bispectrum, 

1/2 

we first note that V" contributes a factor of f^^l ~ (V'" / H) to /nl- We also note that 
what we observe is not f^^°i^, but its projections onto the curvature perturbation (Fig. Sj). 
The efficiency of this projection is determined by the turning angular velocity 6/H. For 
bispectra, we need three transfer vertices (Fig. [3](c)), so the bispectrum is of order 

fNL--P^'^'{e/Hy{V"'/H) . (1.2) 

In this section, we will derive Eq. (11.21) rigorously with explicit numerical coefficients. We 
will also compute the full shape functions using the mixed form of the in-in formalism. 

In Sec. [5l we investigate the shapes of bispectra analytically by taking the squeezed limit, 
Ps Pi = P2, where p^'s are the momenta in the three-point function. We find that the 
three-point function (up to a momentum conservation delta function) scales as 

(Oocp3^/^-% (1.3) 

where u = ■\/9/4 — m? /W^, and m is the effective mass of the isocurvaton. Near = 0, the 
shape changes smoothly to 

(C') (X ln(p3M)P3 • (1-4) 

This scaling behavior lies between that of the local shape ((C^) p^^) and the equilateral 
shape ((C^) oc p^^)- We call these shapes the intermediate shape. We shall explain the 
underlying mechanisms that determine these different shapes. 

In Sec. El we approximate the shape functions by simple analytical expressions, to facil- 
itate future data analyses. 
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Figure 4: The physical interpretation of non-Gaussianity in quasi-single field infiation. The 
non-Gaussianity in the isocurvature (a) direction(s) is not suppressed by slow roll parame- 
ters. This large non-Gaussianity is projected onto the infiationary [6) direction if the infiaton 
trajectory turns. 
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Figure 5: Feynman diagrams for the scalar-exchange four-point correlation function (a) and 
contact- interaction four-point correlation function (b). 
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In Sec. [TJ we discuss the trispectra. As illustrated in Fig. El we show that the trispectra 
are of order 

tNL ~ max |p^-i (e/Hy {V"'/Hf , P^-' (^/i/) V'"j . (1.5) 

Note that ^atl ^ Inl ioi 9/H <^ 1. So in future experiments where angular multipoles / can 
be measured as large as thousands, trispectra may become a better probe for quasi-single 
field inflation. 

We conclude in Sec. [8] and discuss many future directions. 

In Appendix El we present all terms in the Lagrangian up to the cubic order in two 
different gauges, and justify the calculation in the main text. In Appendix |Bl we present 
all terms used in our calculations in terms of the factorized, commutator, and mixed form, 
respectively. In Appendix [Cl we explain the Wick rotation and its several applications in 
this paper. In Appendix |Dl we give general expressions for the series expansion in the in-in 
formalism. 



2 A model 

In this section, we set up an explicit model of quasi-single field inflation evolving two scalar 
fields. The model configuration is illustrated in Fig. [H To consider non-trivial isocurvature 
perturbation while keep things simple, we assume that the parameters characterizing the 
turning trajectory, such as the turning radius and various couplings, are constant. We call 
this the "constant turn" case. 



2.1 Lagrangian, mode functions and transfer vertex 

When the infiaton trajectory moves along an arc, the action is most conveniently written in 
the polar coordinates in terms of the tangential (6) and radial (a) directions of a circle with 
radius R, 



Sn 



(2.1) 



where Kr(6') is a usual slow-roll potential, and Via) is a potential that traps the isocurvaton 
at ctq. The signature of the metric is (— 1,1,1,1). 

The equations of motion for the spatial homogeneous background are similar to those for 
single field inflation. The Hubble equation and the continuity equation are 



R^el ^V,, 



IMlE = R^el . 



(2.2) 
(2.3) 
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The equations of motion for cro(^) and 6o{t) are 

ao = const. , V'{ao) = ROl , (2.4) 

R^Oo + SR^HOo + Kr = . (2.5) 

We have absorbed the constant shift o"o into the net turning radius by defining R = R + ctq. 
The inflaton R6 follows the usual equation of motion (12. 5p determined by the slow-roll 
potential Vsr{0). 

Equation ( 12. 4p indicates that, near ctq, the potential has to contain a linear term to cancel 
the centrifugal force caused by the turning angular velocity. We therefore expand V{a) as 

V{a) = V'{a,){a - a,) + \v"{a,){a - + \v"\a,){a - a^f + • • • (2.6) 

2 D 

with V"{(Jq) > 0{H'^). Depending on the turning angular velocity, the minimum ao of the 
effective potential, 

Kfr = -|(i? + ^-^o)' + na) , (2.7) 

can be far away from the minimum of the original potential V{a). If their separation exceeds 
the typical field range over which the potential can be expanded perturbatively, it is necessary 
to expand the potential around its effective minimum, as we did here. We emphasize that 
the mass V"{ao) and coupling V"'{ao) are evaluated at the new effective minimum ao, and 
are generally independent of their values at a = 0. For simplicity, in the rest of the paper, 
we will just write V, V" and V" and omit their argument a^. 

To study the leading order perturbations, we first present a simple intuitive approach [2] 
and then justify it using a rigorous method. 

We choose the spatially flat gauge in which the scale factor of the metric is homogeneous, 

hij = a\t)5,j , (2.8) 

and the inflaton and isocurvaton are perturbed, 

^(x,t) = ^o(i) + '^^(x,t) , a(x,t) =ao + 5(T(x,t) . (2.9) 

We expand only the matter part of the Lagrangian, namely (12. ip . while ignoring the pertur- 
bations in the gravity part. We neglect terms that are suppressed by the slow- roll parameters, 

_ H _ R'ei Ml ( VL 



iJ2 2if2M2 2 \RVs, 

f V" f V 

— ^ -2M„2-^ + 2Ml — ^ ) . (2.10) 

He ^ RW^r ^ \ RV^r ' 



8 



For cubic terms in the perturbative expansion, we only consider the self- interaction in V{a). 

This approach is intuitively correct because we expect that the main isocurvature contri- 
bution comes from the potential V through the trajectory turning, and in particular, that the 
interactions involving the infiaton field have smaller contribution to the non-Gaussianities. 
Nonetheless, it is important to check this more rigorously, as we shall do in Sec. 12. 2[ 

With the above prescription, we can easily get 

A = - ^im' + pa' - ^id,5aY - pV" - el)5a' , (2.11) 

6C2 = 2a^Reo 6a6e , (2.12) 

£3 = --aW . (2.13) 
6 

£2 describes two free fields, one massless and another massive, in the inflationary background. 
5C2 is the coupling between them, and is the origin of the transfer vertex in Fig. |3] (a). £3 
is the leading source for the three-point function for 66. 
In order to treat SC2 and £3 as perturbations, we need 

This is because the correction to the leading power spectrum is suppressed by the factor 
{60/ H)"^, as we will show; and the requirements that the quadratic term dominates over the 

1 to 



G 



If 



cubic term for 5a < H gives the restriction \V"'\/H <^ (V"/H) . However, it is wort 
mention that these conditions are not the necessary conditions for the model-building 
they are not satisfied, it remains an interesting open question how this model can be solved 
non-perturbatively. 

We define the conjugate momentum densities Siii = dSC/ di^Scpi), where Scpi {i = 1,2) 
stand for 69 and 6a respectively, and 6C stands for the sum of (12.111) . (I2.12p and (I2.13p . We 
then work out the Hamiltonian density in terms of 6(f)i and 611^, and separate them into the 
kinematic part T-Lq and the interaction part Tii. To use the in-in formalism [6], we replace 
50j's and ^vrj's in the Hamiltonian density with the ones in the interaction pictures, 50f 's and 
^tt/'s. These latter variables satisfy the equations of motion followed from the Hq. Finally, 
we replaced 6n^ with 50[ using the relation 6(f)l = &Ho/d{6'7i-). Following this procedure, 
we get the following kinematic Hamiltonian density 

lR'6e/ + $^{dMi? + + ^,{dMif + ImHa] 

2a^ 2 2a'^ 



(2.15) 



^In the absence of higher order derivative terms such as V"" , the \V"'\/H cannot be much larger than 
{V" / HY because otherwise the cr-field approaches a classical instability. But more generally even such kind 
of instability can be rescued by including large higher order derivative terms of V . 
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and interaction Hamiltonian density 



U{ = -C2a^6aiSei , (2.16) 
ni = c^aHa] , (2.17) 



where 



C2 = 2R9o , C3 = IV" , = V" + lOl . (2.18) 

For the constant case, C2, C3 and iv? are all constants, m is the effective mass of the 
isocurvaton|§ 

In the interaction picture, we quantize the Fourier components 56*^ and (5cr^ of the free 
fields 59^ and 5cr^, 

= MkO-k + wlk^-k ' (2-19) 
5<yl = t;k&k + ^* k&lk , (2.20) 

where Ok and 6k are independent of each other, and each satisfies the usual commutation 
relation, 

K, alk'] = {2Tif6\\i + k') , [6k, 6lk'] = {2Tif6\\i + k') . (2.21) 

The mode functions Wk and fk satisfy the linear equations of motion followed from the 
kinematic Hamiltonian, 

2 

< - -< + ^^«k = , (2.22) 

r 

2 

< - + A;V + TT^^^k = 0, (2.23) 

where r is the conformal time, dt = arfr, and the prime denotes the derivative with respect 
to r. 

The solution for the mode functions are 

Mk = —^(1 + ikT)e~'^^ , (2.24) 



and 



= -ie'^-+\)-2^H{-Tf'^Hl^\-kT) , for m'/H^ < 9/4 , (2.25) 



^Although we mostly consider the case m? ^ V" ^ 0{H^) and [Oq/HY ^ 1, we comment that the non- 
perturbative case (Oq/H)'^ > 1 with ^ 0{H^) is also possible because V" can be negative. As long as 
> this does not introduce a tachyon mode at least at the quadratic level. Higher order non-perturbative 
corrections remain an open question. 
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where v = a/9/4 — m'^/H'^; 



= -ie-^2^+^^.^H{-Tf'^H[l\-kT) , for m'/H^ > 9/4 



(2.26) 



where 9 = y^m^/if^ — 9/4. The normahzation of both mode functions have been chosen so 
that, when the momentum k/a is much larger than the Hubble parameter H and the mass 
m, we recover the Bunch-Davies vacuum 



Ru], , vir — )■ I , re 

V2k 



-ikr 



(2.27) 



The behavior of the mode functions after the horizon exit, kr — )■ 0, is also useful. For 



21^-1 



V(-+|)f f_r(z/)^(-r)-'^+i , < z/ < 3/2 



TT k 



1/ = ; 



(2.28) 



for m^/H^ > 9/4, 



fk — )■ — ^e" 



2 



V^l^, .3/2 



(2.29) 



In the kr limit, Eq. ( 121281) contains a decay factor (-r)"''+^/^ while Eq. ( K29\\ has 
both the decay factor (— r)^/^ and an oscillation factor r^*'^. The decay factors indicate that 
the perturbations for massive fields eventually roll back to zero. As the mass increases, the 
behavior of the perturbations changes from that of an over-damped oscillator (corresponding 
to m/H < 3/2) to an under-damped oscillator (corresponding to m/H > 3/2). In the 
under-damped case, the contribution of the massive isocurvaton to the curvature correlation 
functions is suppressed by factors of e"™/^, analogous to the Boltzmann suppression, due to 
the oscillation factor in the mode function. The derivation of this Boltzmann suppression is 
given at the end of Appendix O [Note that this suppression factor is not the first factor e"'^'^/^ 
in fl2.29p . which is cancelled by r(ii> + 1).] Because of this suppression, in the remainder of 
this paper, we shall concentrate on the critical-damped and over-damped case, < u < 3/2 
(corresponding to 3H/2 > m > 0). 



2.2 Two gauges 

In this subsection and Appendix |Al we use a rigorous approach to expand the perturbations. 
This will also justify the simple approach used above. 
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We first still work in the spatially fiat gauge (12 .Sp and f l2.9p . Now we keep all terms, 
including those from the gravity sector, and get 



V" 



3 3 



+ -6a - -{dMY - y (V" - 0'o)Sa' (2.30) 
6C2 =2a^Reoma - 2ea^ReoH6e6a , (2.31) 

5^3 = - —V'"6a' + ■■■ . (2.32) 

6 

The derivation and the full cubic action can be found in Appendix lA.li There we also show 
that the V" term in Eq. f l2.32p indeed gives the leading order contribution to bispectrum. 
Compared with (12. lip and (12.120 . extra terms that are suppressed by slow- roll parameters 
arise in the following two places. Firstly, the 66 field has a small mass, so it is no longer 
constant after the horizon exit. Secondly, there is an extra coupling term between 66 and 
6a, i.e. the second term in (I2.3ip . This term is relatively smaller than (I2.12p before and 
near the horizon exit due to the slow-roll suppression, but can become important afterwards 
because 66 decays exponentially in real time. The detailed form of the coupling is also very 
important to the infrared behavior in the computation of the correlation functions. So we 
would like to clear up these issues by using a different gauge. 

In single field inflation, we know that the quantity 66 is not the conserved quantity after 
the horizon exit, and there are correction terms suppressed by slow-roll parameters when we 
relate it to the conserved curvature perturbation C. This conserved quantity is most explicit 
in the uniform inflaton gauge. Therefore we would also like to use a similar gauge in this 
model and see in particular how the corrections affect the second term in (I2.3ip . We choose 

^(x, t) = 6o{t) , a(x, t) = ao + 6a{x, t) , (2.33) 

and the spatial metric 

hij{^,t) = a2(t)e2'^(^'*)5ij . (2.34) 

Here we have done a position dependent time-shift with respect to the spatially flat gauge, 
so that the inflaton is homogenous but the scale factor has fluctuations. In this gauge, the 
full quadratic Lagrangian is 

A =ea'C' - eaidCr + pa' - \ {d^f - y {v" - dl) 6a' , (2.35) 

6C2 = - 2a^^6aC ■ (2.36) 
H 
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The cubic action is 

£3 = -^y'V + ... . (2.37) 

The derivation can be found in Appendix IA.2I There we also show that the single term 
listed in f l2.37p is the most important one. 

In the uniform inflaton gauge, the physics is conceptually clearer. The linear equation 
of motion for the curvature mode ( describes an exactly massless scalar field in inflationary 
background and has the following leading order solution, 

C, = --^(l + ^kT)e-^^- . (2.38) 

Before considering the coupling to the isocurvature modes, (k is constant after horizon exit 
with H and e evaluated at = aH. The linear equation of motion for the isocurvature mode 
6a describes a massive scalar field, and has the solution f l2.25p and fl2.26p . 6a eventually 
decays away. The transfer vertex is described by the single term in f l2.36p . In our setup, the 
curvature mode ( seeds the large scale structures, as in the single field inflation, while the 
isocurvature perturbation 6a is negligible at the reheating. 

This also gives a justification of the simple method used previously. Comparing (12.350 . 
flOHD and (12381) with fl23ID . (127721) and ([221 , we see that, as long as the parameters H 
and e are approximately constant, we can use the usual time-delay relation 

C ^ ^ -^66 (2.39) 

to convert the correlation functions of 66 to those of C- If H and e are not always constants, 
we should start with ([23SD-(I23ZD. 

We would also like to emphasize that, in our case, although ( is eventually a constant, it 
may not be so right after the horizon exit, in contrast to the single field inflation. As we will 
demonstrate, depending on the mass of the isocurvaton, the isocurvature mode can decay 
slowly and be transferred into the curvature mode long after the horizon exit. 



3 Power spectrum 

The turning trajectory leads to a correction to the power spectrum. The term responsible 
for the transition from the isocurvature mode to the curvature mode is (I2.16p . 

= I rf^x ni = -c,a' I -^M^(^-^ ■ (3-1) 
This term introduces the transfer vertex (Fig. [3]^a)). 
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The two-point function for 59 is given by the expectation value of 59'^ at the end of 
inflation. In the in-in formalism this is given by, 



{59') ^ (0| 



Texp i / dt'Hi{t') 



to 



59]{t) 



Texp 



dt'Hj{t'] 



to 



+ 



to 



dil / dti{Q\Hi{ii) 59] Hi{ti)\0) 



- 2 Re 
+ ■■■ . 



to 



dti / dt2{0\59] Hr{h)Hj{t2)\0) 



to 



|0) (3.2) 

(3.3) 
(3.4) 



The terms (13.31) and (13.41) are the leading order corrections to the two-point function. Ac- 
cording to the Feynman diagram Fig. [3](b), each Hj should be replaced by 

We evaluate these terms using the technique of normal ordering. After normal ordering, 
only terms with all fields contracted survive. For example, a contraction between 5a {p, t) on 
the left with 5a{p',t') on the right gives Vp(t)v*,(t'){27iY5{p + p'). We sum over all terms 
that represent the Feynman diagram Fig. [3](b). 

The term (13.31) gives 



(3.5) 



24u* {t)up,{t) dti a\ti)Vp,{h)up,{ti) (27r)^5^(pi + ps) 



to 



The term (13.41) gives 



Uto 



X (27r)^5^(pi + p2) 



" \Re 



X {2rrf5\p, + p2 



dx x-^'^Hl^'^ (x)e'" (27r)=^5=^(pi + p2 



dtia\ti)vp,{ti)ul^{ti) I dt2a\t2)vl^{t2)ul^{t2) 

J to 



dxi x^^^^Hl^\xi)e-'''' / dx2 x^^'^H'^\x2)e 



(3.6) 



For 1/2 < V < 3/2, there are divergences at r — t- in Eqs. (13.51) and (13.61) . However 
these divergences are spurious and are canceled by summing up (13.51) and (13. 6p . The final 
result can be written as 



(2vr)V(pi + p2) 



,3 ' 



(3.7) 



where the numerical factor C is defined as 



Civ) = -Re 
4 



/ dxr / dx2 [x^^'^H^J\xi)e''''x:,^'^H'^\x2)e-'''^ 

Jo Jxi ^ 



x;^^^Hl^\xi)e-'^'x:,'''Hi;'^{x2)e 



ix.-l/2TT{2), 



-1X2 



(3.8) 
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Figure 6: The coefficient defined in Eq. fl3.8p . 

It is not difficult to see that, as X2 — )■ Xi — )■ 0, the leading divergence in the integrand is 
imaginary and C(z/) is finite for u < 3/2. We plotted C as a function of u in Fig. [6l 
We use fICTD to relate (C^) to (56^), and get 



(0 = (27r)V(pi + p2)-3Pc, 



where the power spectrum is 



Pr 



The spectral index is 



(3.9) 



(3.10) 



1 



dlnPr 



C 



dink 



-2e-7] + SCt] I 



(3.11) 



The last terms proportional to {Oq/H)"^ in f l3.10p and (13. lip are the corrections from turning 
trajectory. So the perturbation theory requires {Oq/HY <^ 1. 

If the turning of trajectory only happens for a period of time, transient large running of 
the power spectrum and spectral index becomes possible. We leave this non-constant turn 
case to future work. 

Although the observational consequence on the power spectrum for the constant turn case 
is not very distinctive, this case serves as a warmup exercise to the bispectra and trispectra 
cases in the next sections. As we will see, the simple versions of the problems and solutions 
that we have encountered here will be significantly magnified. 
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4 Bispectra 



The origin of the large non-Gaussianities is the self-interaction terms for the isocurvaton 
field a. The leading cubic term in this model is fl2.17p . This gives the following term in the 
interaction Hamiltonian, 



(4.1) 



Through the transfer vertex (Fig. [3]^a)), this interaction is converted to that of the curvature 
perturbation. The Feynman diagram is given in Fig. El^c). 
The three-point function for 66 is given by, 



(0| 



Texpii / dt'Hi{t' 



to 



69Kt) 



Texp 



dt'Hiit'] 



to 



|0). (4.2) 



Again using fl2.39p . {(') is related to it by a factor of {—H/6q)'. 

One can expand (14. 2 p and write it in two equivalent forms. In the first form, we simply 
expand the exponentials. 



{69') 



to 



to 



to 



dil / dk / dti / dt2 {Hi{i2)Hi{ii) 69] Hi{ti)Hi{t2)) 



to 



- 2 Re 
+ 2 Re 



to 



to 



to 



t2 



dh / dh / dh / dh {Hj{h) 69-j Hj{h)Hj{t2)Hi{ts)) 



to 



Uto 



to 



t2 



to 



ts 



dh / dh / dt3 / dU {69j Hj{U)Hj{t2)Hj{t,)Hi{U)) 



to 



(4.3) 
(4.4) 
• (4.5) 



In each term, the interaction vertex can appear in one of the four Hj, and the rest three 
i?/'s should be replaced by the transfer vertex We refer to this form as the factorized 
form because there is no cross time-ordering between the integral from the left interaction 
vacuum and that from the right. After contractions, we get 



X Re 

X 
X 



Tl 



dri a {ti)v* {ti)u (ti) / dT2 a {T2)Vp^{T2)Vp2{r2)Vp:^{T2) 



dn a\Ti)v* {Ti)u*{ri) 



n 



dT2 a {t2)v*p.^{t2)u*^{t2] 



(27r)^5^(^^ Pi) + 9 other similar terms 

i 

+ 5 permutations of pi . 
We leave the details of the 9 other terms to Appendix IB.li 



(4.6) 
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In the second form, ( I4.2p can be expressed in terms of the nested commutators [H], 
{66'')= fdt, r dh r dh r dU{[Hj{U),[Hj{t3),[Hj{t2),[Hj{ti),6ej{t)']]]]) , (4.7) 

J to J to J to J to 

in which all the integrands are written under a single time-ordered integral. We refer to this 
form as the commutator form. Again replacing one of the Hj with and the rest with H2, 
we get 

3) = UclcsUpMupMupM 



X Re 



/O PTl PT2 pT3 J* 

dn / dT2 / dT3 / dTi Y\ 
-00 J —00 J —00 J —00 •_1 



X {a{T2)A + a{r,)B + a{r^)C)] {27rf6'{J2Pi) 

+ 5 perm. . (4-8) 

We have used the fact that Up. {0) are real. A, B and C are three contributions corresponding 
to replacing Hj(t2), Hjit^) and Hjit^), respectively, with 

^ = Ki(^l) - C-^-) (^Pi(^l)<(^2) - C.C.) {Vp.^{T2)vl.^{Ti)u'*^{r^) - C.C.) 

^P.(r2)t;;(r3)<(r3) , (4.9) 

^ = Ki(^l) ~ C-C-) K2(^2) - C.C.) {vl^{Ti)vl^{T2)Vp^{T-i)Vp^{T:>,) - C.C.) 

^P3(^3X(r4X(r4) , (4.10) 
C = - «(ri) - C.C.) K^(r2) - C.C.) («;3(r3) - c.c.) 

<(^iX(^2)f;3(r3)t;pi(r4)fp2(r4)fp3(r4) . (4.11) 

The factorized and commutator form each has its computational advantages and disad- 
vantages when evaluating the integrals. To see this, let us investigate the properties of these 
integrals in the IR (r — )• 0) and UV (r — > —00). 

• IR convergence. To see the IR behavior, we use the asymptotic forms of the mode 
functions in the r — t- limit, 

u'pM) oc (-r)(l-ipir + ---) , 

VpXr) (X {-r)l-'^ {I + a,{-rf + a2{-rr + ■ ■ ■) , (4.12) 

where we have ignored an overall phase in Vp^ that will always be cancelled. Also note that 
ai is real, ^2 is complex. 

We first look at the factorized form. It is easy to see that each term in (14. 6 p has IR 
divergence for 3/2>z/>l/2(0<m< \/2H). For smaller m, the isocurvature mode is 
decaying slower, hence the conversion to the curvature mode lasts longer after the horizon 
exit. But since the isocurvature mode eventually decays, we do not expect any singular 
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behavior starting from m = \/2H, as we have seen in the case of the power spectrum. So 
these IR divergences should be cancelled when we sum over all ten terms. It is possible to 
check this analytically, but it takes very long time even with the aid of Mathematica. This is 
because we not only have the leading order spurious divergence, ~ t^'^", but also have eight 
different orders of subleading spurious divergence, such as ~ t^'^" , t^~^^ , ■ ■ ■ . Numerically, 
errors occurred in the cancellation of these spurious divergence quickly dominate in the final 
results as z/ approaches 1/2 from below. 

However this cancellation is made much more transparent in the commutator form. In 
this form, all integrands are under the same integral and their mutual cancellation in the 
IR is explicit. Without loss of generality, let us examine the term A fl4.9p . Using fl4.12p . 
we see that if we take Tj — )• 0, the leading powers of in the whole integral go as follows: 
(-ri)2/2^^, (-r2)^/^~^'', (-rs)^/^-!^^ (^_^^y/2-v _ Yqi the multi-layer integral (gT]), if an inner 
integral diverges in the IR, its contribution to the outer integral is dominated by this IR 
behavior; if it converges in the IR, its contribution to the outer integral is 0{1) and complex. 
One can then list all possibilities and examine them case by case. But a quicker way to see 
the conclusion goes as follows. The largest power for Xj that we listed above is 3/2 — z/. So 
by considering the case v = 3/2, we are considering the largest possible IR contributions 
for all < < 3/2. For this case, we can just take the Tj — )■ limit for all terms in the 
integral and take the limit Tj — )■ ti. The terms in the first bracket in fl4.9p goes as (— ti)^ as 
the leading term in is cancelled by the subtraction of its complex conjugate. The terms 
in the second and third brackets are similar, they go as {—ti)^ and (— ri)^"^*^, respectively. 
Since the first line of (14.91) is pure imaginary, the second line, Vp^v*^u'*^, has to be imaginary 
to make the overall integrand real. Hence it goes as (— ri)^"^'". Including all the scale factors, 
a{Ti) (X 1/Ti, the whole integral goes as (— ti)^~^'^ = 1 and is logarithmically divergent for 
u = 3/2. So we conclude that the IR divergence is explicitly absent for u < 3/2 in the 
commutator form. 

• UV convergence. It is a common feature that the integrands in these correlation func- 
tions are oscillatory in the UV, when modes are well within the horizon. Their contribution 
is averaged out. For the Bunch-Davies vacuum, this regulation can be achieved by slightly 
tilting the integration contour into the imaginary plane, r — > r(l ± ie). 

For the factorized form, we tilt the contour clockwise, r — > r(l — ie), for the anti-time- 
ordered integral from the left, and counter-clockwise, r — )■ r(l + ie), for the right. This 
procedure works for all u. 

For the commutator form, if we do not have the problem of the spurious IR divergence, a 
procedure similar to the above still works. Now the left and right factors are mixed and the 
original tilts cannot be kept intact, but the effect of these tilt can be easily implemented. The 
effect is to suppress the oscillating contributions in UV. So we can re-choose a proper tilt for 
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view of the factorized form 

i 



each term in the new integral according to its convergent direction. The only exception occurs 
for some special momentum configurations where the oscillating factor form the left and 
right happen to cancel each other. At these places, one encounters the slower convergence, 
or spurious divergence that are completely absent in the point o 
discussed above. These divergences eventually will be cancelled 

However with the problem of the spurious IR divergence, as for the case l/2<z/<3/2, 
the problem in UV gets much worse. We cannot even avoid the spurious UV divergence 
for general momentum configurations, while preserving the explicit IR convergence. This 
is because, for example after expanding the term A (14. 9p . different terms have different 
convergent tilting directions as we mentioned. But now we cannot re-choose the tilts for 
them individually if they have to be grouped to achieve the explicit IR convergence. 

• Mixed form. In summary, the factorized form is much more convenient to achieve the 
explicit UV convergence, while the commutator form is much more convenient to achieve the 
explicit IR convergence. For < z/ < 1/2, we do not have the IR problem, so we can use the 
factorized form, but it starts to fail as z/ — )■ 1/2. For 1/2 < u < 3/2, both types of spurious 
divergence exist. The best way to proceed is to combine the two forms. We introduce a 
cutoff Tc (e.g. Tc = —2/pi), and write the IR part (r^ < r < 0) of the integrals in terms of 
the commutator form, and the UV part (r < Tc) in terms of the factorized form. 



rn-1 ptc rTn-i 

/ dri ■ ■ ■ I dr-i {commutator form} / dTi^i ■ ■ ■ dr^ {factorized form}(4.13) 

^ J Tc J Tc J —CO J —CO 

The detailed expressions of this mixed form are presented in Appendix IB.3[ 

• Wick rotation. An additional subtlety is that the contour tilting is easy to do analyti- 
cally, but difficult to implement numerically. In this paper, we will use a much more efficient 
approach of Wick rotation to achieve the fast convergence in UV. As shown in Appendix 
\C\ one can rotate the integration in the complex plane, Tj — t- ±ixi, so that the oscillation 
factors for — > — oo become the exponential suppression factors in the Xj-coordinates. 

Applying the technique of Wick rotation to the mixed form, we finally have a very efficient 
way to numerically compute the full shapes of the bispectra. As an example, we present the 
plots for 1/ = 0, 0.3, 0.5, 1 in Fig. [3 

•^To illustrate this using a simple example, we look at the integral dxie~'^°'^^ ■ j'^^dx2e^''^^ = 
J^^dxi (ia;2(e~"'^ie*^^^ + e~*''^2g«&a;i^_ Tilting the integration contours at x — ?► — oo into the imag- 
inary plane according to the signs of a, b and a ~ b, the factorized form on the LHS simply gives l/a&, and 
the commutator form on the RHS gives l/(a — 6)6 — l/(a — b)a. Each of the two terms in the commutator 
form has a divergence at the special point a = 6, but cancelled by each other. Cases of such UV spurious 
divergences are encountered, for example, in the study of trispectra [SUllj if one uses the commutator form. 
So the factorized form is more convenient in such cases. 
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Figure 7: Shapes of bispectra with intermediate forms. We plot {piP2PzY F with v = 
0, 0.3, 0.5, 1. The plot is normalized such that {piP^V^YF = 1 for Pi = P2 = Ps = 1- 

To plot, we define the function F as 

= F(pi,P2,P3)Pc'(2^)'^'(Ep^) • (4.14) 

i 

To illustrate the shape of a scale-invariant bispectrum, we conventionally normalize the 
amplitude F by multiplying a factor of {piP2PzY- This makes it dimensionless and scale- 
independent. 

From Fig. [71 we can see that when v is small, the shape looks more like an equilateral 
shape. When v gets larger, the shape looks more like a local shape. In Sec. [5l we will study 
the analytical properties and explain the underlying physics of these shapes. 

Finally, we would like to parameterize the magnitude of the non-Gaussianities in terms 
of an estimator f^^\. According to the convention in the bispectrum literature, we define the 
number by matching with the /j^^*^^ in the local shape ansatz in the equilateral limit. In 



20 



a(v), 



300 


: 2.0; 


250 


: 1.5- 


200 


r 1.0 - 


150 


: 0.5 1 


100 





50 





0.0 



0.2 



0.0 0.2 0.4 0.6 



0.8 



0.4 



0.6 



o.e 



1.0 



1.2 



Figure 8: The numerical coefficient a(z^) in 



int 
NL- 



other words, we define /]^^ so that, at pi = P2 = Ps, the three point correlation function is 



(C(Pl)C(P2)C(P3)) ^ (2vr) V(pi + P2 + P3)Pc' ( Y^fNL 



9 



int 



1 

7i 



From this definition, we get 



^int 

Jnl 



a{u)p-"^[do/H] {-V"'/H) 



(4.15) 



(4.16) 



verifying the qualitative estimate given in the Introduction. The numerical coefficient ot{v) 
is plotted in Fig. [HI and ^ 6.1 x 10~^. Oq is positive in our convention, so the sign of the 
is the opposite of the V"'{aQ). 
From Fig. [HI we can see that «(z/) can get very large. For example, as v varies from 
to 1.35 (i.e. mass varies from l.bH to 0.65-/^), a(z^) grows from 0.2 to 300. As z/ — 3/2, 
a{v) blows up. This divergence happens because we are using the constant turn assumption. 
When the effective mass of a becomes zero, a fluctuation 5a never decays at super-horizon. 
Then the transfer from isocurvaton to curvaton lasts forever. Practically, the upper bound of 
the conformal time integration will not be zero. If the horizon crossing time of a perturbation 
mode is Nf e-folds before the end of inflation (or the time when the inflaton trajectory 
becomes straight), one needs to impose a cutoff 

As 1/ — )■ 3/2, one can show that (14.81) is dominated by the integrals that behavior as 
ni=i / dTi/Ti ~ Nj. In principle, Nf can be as large as 60. But we caution that, al- 
though a{v) grows a lot as — )• 3/2, this does not mean that the non-Gaussianities can 
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be enhanced by such a large factor, at least in the perturbative regime. The reason is the 
following. In this limit, in (13. 8p scales as Nj for the same reason. For large Nf, the 
perturbation theory requires Nj{9Q/ H)"^ <^ 1 instead. Therefore, in the perturbative regime, 
the effective enhancement factor is only Nf. 



5 Squeezed limit of bispectra 

In this section, we investigate the behavior of the three-point correlation function in the 
squeezed limit, ps <^ pi = P2- This is important for several reasons. First, in this special 
limit, analytical results for the shape functions are possible. These provide both useful 
checks on our numerical results and complimentary information. Second, numerical results 
are not useful in the construction of estimators that are often used in data analyses. We 
need to guess simple analytical expressions if they are not immediately available. Knowing 
the analytical results in the squeezed limit greatly helps in achieving this goal. Third, the 
scaling behavior of the squeezed limit is closely tied to the underlying physical mechanisms 
and we will use it to classify the shapes of bispectra. 

We will investigate this limit in both the commutator and factorized form. We will see 
that they give equivalent results. The following behavior of the Hankel function in the small 
argument limit, x ^ 1, will be useful in the analyses, 

, 2^(1^) „ 2-2+^r(i/) / cos(7ri^)r(-z/) 1 \ , , 

hI^\x) ^ -I ^-^x"" - i , , [ x-"^^ + -I ^-^^ + —— x'^ + ■ ■ ■ (5.1) 



Note that the real part starts from 0{xf\ 

We start with the commutator form and first look at the contribution from the A term 

^3 ^3 



SVT'^ C^C3 1 



X Re 



26 HR^pip2Pz 

' dn r dr, r dr, r dn {-r{)-'/\-r,f'\-r,)-''\-r,r''' 



X 
X 



sin(-piri) [H^^){-p,r,)H^;^\-p,T2) - c.c.) {H^;^\-p^T2)H^J\-p^T,)e-'P-^^' - c.c.) 
H^^\-P2T2)H^^\-P2n)e'^'^'] + 5 perm. . (5.2) 



We neglect the common factor (27r)^(5^(^ pj) in this section. Terms with momentum permu- 
tation behave differently in the squeezed limit, and we exam each of them in the following. 
For the term exphcitly written in (15.21) . we define = piTi {i = 1, 2, 3, 4) and get 

Stt^ c|c3 1 



26 HR^pfp2P3 
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X Re 



dxi 



Xl 



^2 



^3 



dxd 



-x,)-"\-x,fl\-x,)-''\-xA-''' 



sin(-a;i) [hI'){-x,)H':^\-X2) - c.c.) ( Hi'\-^X2)Hi'\-^x,)e 

\ Pi Pi 

Hi'\-^-^x,)Hl^\-P^x,)e''^^^^ 
Pi Pi 



PI — C.C. 



(5.3) 



The terms H^\—X2P'i/pi) in the 3rd hne can be approximated in the small —x^Pzjpi limit. 
The reason is as follows. If —X2Pi/pi ~ 1, \x2\ ^ 1. Other terms in the integrand have factors 
such as H^^\x2). These factors become fast-oscillating and hence suppress the integration. 
However the terms H^\—X4,p^/pi) and e"*^*^^/^^ in the 3rd line cannot be approximated in 
the small —Xi^p^lpi limit, since there is no oscillatory term as X4 gets large. We redefine 
?/4 = Xip^/pi. With this prescription we get 



26- 



-r z/ 







X Re 



Pi P2 Pi 



dxi 



Xl 



dxi 



X2 



dx-^ 



PI 



X3 



dyA 



-xi)-'/\-X2f'\-x,)-'/\-y,)-'l^ 



X sin(-xi) {H^^){-x,)H^^\-X2) - c.c) {-X2)-' H^^\-X2)H^^\~x,y 
X {H^^){-y,)e-^y^ + c.c.)] . 



(5.4) 



In order not to be suppressed, 0:3 < 1. So the upper limit of the integral is effectively 0. 
Since it is also convergent at ?/4 — > 0, this integral can be factored out. 

We next look at the term with the permutation pi Ps. With the same definition of Xi 

{i = 1, 2, 3, 4), we get 

37r3 c^c3 1 



26 HR^pfp2P3 



X Re 



dxi 



Xl 



dx-) 



X2 



dx-i 



^3 



dxA 



-X,)-"\-X2Y'\-X,)-"\-X,)-'/' 



il)(_P3 
Pi 



X sm( Xl 

Pi 



Xi)Hl'\-^^X2) - c.c. j {Hi'\-X2)Hi'\-x,)e-''^ - c.c.) 



X X2lir{ XsCn 

Pi Pi 



(5.5) 



In the third line, the first three functions can be approximated in the small —Xip^/pi {i = 1, 2) 
limit, and one of the Hankel functions should be expanded to 0{x'() in order to get a non-zero 
result. Focusing on the scaling behavior of pi, we see that it is proportional to 

1 



P1P2 



(5.6) 



Comparing to (15. 4p . this is negligible. 
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Finally, we look at the term with the permutation p2 -H- pa, 
Svr^ clcs 1 

dx2 / dXs / dX4 {-Xi)-^/\-X2y/\-X3)-^^\-X4y^^^ 
oo J —oo J —oo 



X sm 



V Pi Pi 



Pi Pi 



(5.7) 



In this case, we can also approximated the three functions in the 4th line in the small 
—Xip^/pi {i = 2, 3) limit. For u > 1/2, we use the leading term for the two Hankel functions 
and the subleading term for the exponential function and get 

(5.8) 



pi-"'P2Pr ' 



and for u < 1/2, one of the Hankel functions should be expanded to 0{x^) and we use the 
leading term for the exponential function, 

~ -I^ ■ (5.9) 

PiP2P3 

Both (15. 8 p and (15. 9p are negligible comparing to (15. 4p . for z/ < 3/2. 
The other permutation pi ^ p2 gives each term a factor of 2. 

We perform the similar analyses to the B and C terms. Overall we find that the dominant 
contribution come from the A and B terms and their momentum permutation pi O p2 only. 
The final result is 

(C(Pi)C(P2)C(P3)) ^^^^^ ^(^)|| . / 3^, (2vr)^^^(EP-) ' (5-10) 

Pl P2 Pi i 



where 



37r^r(z/) 

dxi / dx2 / dx3 

oo J —oo J —oo 



[(-Xi)-V2(_a;2)i/2-(-a;3)-^/^in(-Xi) 
Im (Hl'\-x,)Hi'\-X2)) Im {Hi'\-X2)Hl'\-X3)e''-^) 
+ (-Xi)-i/2(-X2)-'/'(-X3)'/'-'^ sin(-xi) sin(-X2) 
Im (Hl'\-x,)Hi'\~X2) {Hl'\-X3)y' ' 



X 



r dy,{-y,r'/'Re {Hi'\-y,)e-'y^) . (5.11) 

•J —oo 
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Similar analyses can be applied to the ten terms of the factorized form listed in Appendix 
[RTl We find 



25-u 



X Re 



—I 



X 



+ i 



dxi 



-x,)-'/\-x,y/'-^Hl'\-x,)e-^ {Hi'\-x,)y 



-X2 



-'/'Hl'\-x,)Hj^\-xM'\-x,)e-^) 



dxi {—Xi] 



-1/2^(2) (_ 



-Xi)e 



dxi 



oo 




dxo 



-Xi 



1-1/2 



{-x,)-'/'H^J\-x,)e-^H^^\- 



X2 e 



1X2 



dxi 



XI 



dxo 



X2 



dx'^ 



-X2 



-oo 
1-1/2/ 



(2), 



^1X3 



1X2 tt(2) 

V 



ixz 



+ (-x,)-'/'i-X2)-'^'{-XsY^'-^Hi'\-x,)e'''^Hl'\-X2)e''' (hI'^-xs)} 
X / d£2(-X2)"'/'Re (//('H-£2)e-*"^) . (5.12) 

J — oo 

One can rewrite the expression Re[---] in fl5.12p in terms of one time-ordered cubic 
integral, and this exactly reproduces (15. lip . Therefore the two expressions are equivalent. 

The computational advantage of each form is the same as we have discussed in Sec. [3 
The fastest way to evaluate s{v) is to combine them into a mixed form, and do a Wick 
rotation in the UV part of this form. In this special case of the squeezed limit, it turns out 
that the expression (15.1 ip can also be evaluated by brute-force without these treatments. 
This is because the cubic integral is computationally less expensive than the quartic integral, 
and the UV behavior of this cubic integral is roughly J dx x~'^~^/^e^*^. For 1/2 < u < 3/2, it 
converges without any regulation; for < z/ < 1/2 the oscillatory factors help it to converge, 
although the speed is increasingly slow towards = 0. We plot s{v) in Fig. |9l 

The squeezed limit as z/ — )■ needs some extra care. When we use (15. ip in the above 
analyses, the subleading term is suppressed by a factor of {pz/viY^ ■ This is the case only if 



Pi 



(5.13) 



If V is sufficiently close to zero so that the above condition is no longer satisfied, the expansion 
of the Hankel function should be changed from 



Pi 



Pi 

Pi 



(5.14) 
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Figure 9: The coefficient s(z/), defined in Eq. (15 .111) . 



to 



Pi VT Pi 

If we fix ps/pi while reducing z/, as z/ < — (ln(p3/pi))^^, the factor 

37r2r(z/) 1 



Pi P2 Pi 



in f lS.lOp and (15.111) should be changed to 



2^ HR^ A A 
Pi P2 P3 



(5.15) 



(5.16) 



(5.17) 



and the two factors of {—Xi)~'^ inside the integrals in (15. lip should be changed to —1. In 
Fig. ini we have assumed (I5.10p . Therefore, the divergence as z/ approaches does not mean 
that the non-Gaussianity is blowing up, rather signals the change of shape. 

We end this section by discussing two interesting physical aspects. 

• Quasi- equilateral and quasi-local shapes. Notice that the squeezed limit of our bispectra 
(~ Ps^^'^~'^) lies between that of the equilateral (~ p^^) and local (~ pj^) bispectrum p^[T3]. 
We call these shapes the "intermediate shapes" . They are not super-position of any previous 
known shapes. For example the superposition of the local and equilateral shape gives a 
different scaling behavior at the squeezed limit. For the constant turn case, these bispectra 
are scale- invariant. 

Equilateral and local bispectra are two well-known types of scale- invariant non-Gaussianities 
that can become observably large. The underlying physics associated with these two shapes 
are as follows. 
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The large equilateral bispectrum is typically generated when the interacting modes are 
crossing the horizon around the same time. This happens for example in the single field 
infiation models with higher derivative interactions or certain multifield generalizations [13]- 
[T6] . such as DBI inflation [TTtfTS] or k- inflation [T9||20]. Long wavelength modes that already 
crossed the horizon is frozen in single fleld inflation, and cannot have large correlations with 
the modes that are much shorter. As a consequence, the bispectrum, properly normalized, 
peaks at the equilateral limit. 

On the contrary, the large local bispectrum is typically generated when the modes have 
already exited the horizon. This happens for example in special types of multi-fleld slow- 
roll models [211[22] or curvaton models [23]. Superhorizon curvature perturbations are not 
conserved in multi-fleld inflation models, and can receive contributions from isocurvature 
modes. Different patches of universe that are separated by inflationary horizons evolve 
independently, and so the non-Gaussianities come in locally in position space. Thus in 
momentum space the correlation becomes non-local and peaks in the squeezed triangle limit. 

In the quasi-single fleld inflation, the isocurvaton has a mass that can vary around 0{H). 
For the heavier fleld m > \/2H^ i.e. v < 1/2, its amplitude decays faster after horizon-exit. 
Therefore large interactions happen during the horizon exit, and we get shapes that are 
closer to the equilateral type. Namely, the properly normalized bispectra peak roughly at 
the equilateral limit We call it "quasi-equilateral". The numerical example of such a shape 
can be found in Fig. [7] (i^ = 0.3). 

For the lighter fleld m < \/2H, i.e. 3/2 > u > 1/2, its amplitude decays slower. The 
conversion from the isocurvature to curvature mode is still continuing after the mode exits 
the horizon. Here the interactions among the isocurvature modes are local, and in addition 
the conversion will become increasingly local in the position space for the reason that we 
have explained. So we get shapes that are closer to the local type, and they peak at the 
squeezed limit. We call it "quasi-local". See Fig. [7] (z/ = 1). 

In the limiting case = 0, i.e. u = 3/2, the squeezed limit of our bispectra coincide 
with that of the local type. In this limit, the isocurvaton fluctuations do not decay. This 
is the reason that, in Fig. |6l |8] and [9l the amplitudes of C, f^^\ and s approach inflnity as 
z/ — 7- 3/2 for the constant turn case. As discussed in Sec. HI infrared e- folds cutoff should be 
considered in this limiting case. Most of our analyses still apply in this limit, but we would 
like to distinguish the following two cases. First, if V" is still large, we can use f l4.16p but 
with infrared e-folds cutoff. The cutoff will introduce a running in f^^)^ because different 
modes correspond to different Nf and a(i/) is iVj-dependent. Second, if the potential in the 
isocurvature direction also becomes a slow-roll potential, V" is very small, ~ 0{e^^'^)H'^ / Mp. 

''The fact that the shapes are in general flatter than the equilateral shape has to do with the fact that 
the interaction in this model originates from V{(j) and is local in the position space to start with. 
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In this case, other terms in the cubic Lagrangian will become more important, but in terms 
of contributing to /j^^ they are all small. As we discussed in Sec. HI the enhancement 
factor from a(z/) is only Nf in the perturbative regime. So the bispectrum in this case is 
small. This is in accordance with the general findings of previous studies that it is very 
difficult to generate large non-Gaussianities in terms of turning trajectories in multifield 
slow-roll models [22], essentially because imposing slow- roll conditions in all directions are 
very restrictive. 

As we have seen, the shapes of the non-Gaussianities depend very sensitively on the mass 
of the isocurvaton. The shape of bispectrum changes from quasi-equilateral to quasi-local as 
m changes just from lAYH {u = 0.3) to about I.IH {y = 1). However on the other hand, 
as u becomes close to 1.5, the shapes of bispectra are very close to the local form, but m is 
still of order H. For example, for z/ = 1.4, m is still ~ 0.54iJ; but the shape scales as p^'^'^, 
comparing to the local one pj^. It is clear that m is still much too heavy for the quasi-single 
field infiation to become two-field slow-roll inflation, hence the underlying models are very 
different. Therefore, it is an interesting question how good we can distinguish the quasi-local 
form from the local form experimentally. 

• Change of shapes in isocurvature-curvature conversion. The non-Gaussianities in this 
model originate from non-Gaussian fluctuations in the isocurvature direction. It is interesting 
to look at the shapes of the three-point correlation function of the isocurvature modes before 
it is transformed into that of the curvature modes, 

J to 

= ^c■J^H\-rf/'Hl^\-p,r)H^:^){-p,r)H^^^ 

f dn{-ny/'Hl'\-prn)Hl'\-p2n)Hi'\-p,n) + c.c. . (5.18) 

In the squeezed limit and for modes that exit the horizon, this is proportional to 

~ (-pir)9/2-3>3>-6+2. _ (5,19) 

We see that its amplitude is decaying and its shape goes as p^'^'^- Therefore it is evident that 
the effect of the transfer vertex is not a simple projection. During the transfer, the shape of 
the correlation function has been changed, slightly towards the local type. It is important 
to investigate such changes in other cases including the multi-fleld inflationary models. 

6 Shape ansatz 

As we have seen, the precise shapes of the bispectra are complicated for quasi-single fleld 
inflation and we have to perform numerical integration to see the full shapes. However, for 
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the purpose of data analyses, simple analytical expressions which resemble closely to the 
precise shapes are desirable. So we would like to start with the analytical results in the 
squeezed limit, and construct such shape ansatz. For example, we find the following ansatz 
reproduces the squeezed limit behavior and has good overall match with our numerical 
results, 

. - /Jfi (^miEL^ , (6.1) 



where Ny is the Neumann Function and F is defined in (I4.14p . The parameter a can be 
adjusted to fit the ansatz with the numerical results, we found a ~ 8. For data analyses, a 
is no longer a free parameter. 

We have compared this ansatz with the numerical results of the full shapes for v = 
0, 0.2, 0.3, 0.5, 1, and found good match. Examples of the shape ansatz are shown in Fig. [TOl 
and should be compared with the numerical results in Fig. [71 The overall shapes as a function 
of u and the squeezed limits of both results match well. There are some small differences in 
the bulk. For example, for the z/ = 1/2 case, although there is some slight growth towards 
the folded triangle limit in the numerical results, we have numerically checked that the shape 
remains very flat if we look at the more squeezed configuration, consistent with the analytical 
results. So these small differences should be due to the second order terms that we did not 
take into account in the analytical computation of the squeezed limit. 

In data analyses, the construction of the estimator involves triple integral of the shape 
function over the three momenta pi. To have practical computational costs, it is necessary 
to factorize this integral into a multiplication of three integrals that involve individual pi . So 
we would like to further approximate the shape ansatz by templates with simpler functions, 
for example, 

J. 3§-=^- f^iipj+pl+pl) 

10 (piP2P3)2"^''(pi + P2 + Ps)"^'^" 

These simpler templates are shown in Fig. [TTJ They reproduce the shape functions quite 
well except near u = 0. If necessary, one can come up with other templates that are more 
factorized, to get rid of {pi + p2 + Ps)''^^'^'^^'^ ■ 

7 Trispectra 

Four-point correlation functions provide complementary information on inflationary dynam- 
ics, as well as redundant checks on the three-point functions. The study on the trispectra 
has recently attracted much attention in data analyses [21] and model building 
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Figure 10: Shape ansatz ( 16. ip with a = 8. We plot {piP2P3)'^F for i/ = 0, 0.3, 0.5, 
respectively. The plot is normalized such that {piP2'PzYF = 1 for Pi = P2 = Ps = 1- 
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Figure 11: Shape templates (16 .2^ with the same convention as in Fig. [101 
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In our model, there are two terms in the interacting Hamiltonian that contribute to 
the leading trispectra, namely, the V" and V"" terms. There are two Feynman diagrams, 
the contact-interaction diagram (Fig. EU^a)) and the scalar-exchange diagram (Fig. Et^b)). 
Analogous to the fiy^ in the bispectra, we use the t^i to denote the magnitude of the 
trispectra for a given shape. For each shape component, we take the regular tetrahedron 
limit and define t^L as [9] 

(C')component ^ (27r)9p3^3( V p,)\tNL • (7.1) 
limit ^ ^—^ 

i 

For the contact-interaction diagram, the contribution from the V"" term can be calculated 
by expanding the in- in formalism to the fifth order (Appendix [D|) . 



{Se^) D -2Re 

-|-2Re 
-2Re 



dii / dh / dU / dt2 / dh {Hj{i2)Hj{ii) Sej Hj{h)Hj{h)Hi{t3)) 

to •'to ''to J to J to 

rfti / dh / dh / dh / dh {Hiih) Se^ Hi{h)Hj{h)Hj{h)Hj{h)) 

to •'to •'to J to •'to 

ti rt2 pts rti 



dh r dh dh r dh dh {6ej Hi{h)Hi{h)Hi{h)Hi{h)Hi{h)) 

Jto J to Jto J to 

(7.2 



or alternatively using the commutator form 

{59^)^1 f dh-- - r dh{[Hj{h)AHi{h),--- ,[Hi{ti),Qi{t)]---]]) , (7.3) 

Jto •'to 

where one of the Hj's is to be replaced by Hi = ^ J d^xa?V""5a^, and other Hjs are to 
be replaced by The order of magnitude estimate for the size of this trispectrum is as 
follows. The quartic interaction contributes a factor of V"" to tTVL- The transfer efficiency 
is ~ {6/H)^, because we need four transfer- vertices. Keeping in mind that ( ~ V^^'^ in (C^), 
and comparing with the definition (17. ip . we get 

t%\r^P^-'{e/HyV"" , (7.4) 

where the superscript "CI" denotes contact-interaction. 

For the scalar-exchange diagram, the contribution of the V" term can be calculated by 
expanding the in-in formalism to the sixth order (Appendix |D]) , where two of the if/'s are to 
be replaced by H^, and other if/'s are to be replaced by i/g- Each cubic vertex contributes 
a factor of {V"'/H) to t^L, and the transfer efficiency is again {6/H)^. So overall, 

t%\-P^-'{e/Hy{V"'/Hf , (7.5) 
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where the super-script "SE" denote scalar-exchange. 
Comparing f l7.5p with the bispectra fl4.16p . we have 

tfL-{H/eyf^,. (7.6) 

In the slow turn case, {o / ^ 1, so tjvL ^ Inl- Such a large trispectra may be a better 
probe for quasi-single field infiation than bispectra. Comparing (17. 4p with the bispectra 
(I4.16p . we have 

t''^L-{Hie)\v""H'/{y'"f)pNL- (7.7) 

t'^]^ can be either larger or smaller than /^^ and depending on the details of the 
potential V" and V"" . 

It is important to study the integration numerically and analytically to see the shapes of 
the trispectra, as well as the coefficients in front of Eqs. (17. 4p and (17. 5p . This is beyond the 
scope of the current work. 



8 Conclusion and discussion 

To conclude, in this paper, we have investigated in detail a quasi-single field infiation model. 
We find fields with mass of order H can have important impacts on density perturbations 
through, for example, turning trajectories. These effects can be computed perturbatively 
using transfer vertex and Feynman diagrams in the in-in formalism. A one-parameter family 
of potentially observable large bispectra arise. These new shapes are controlled sensitively 
by the mass of the isocurvaton and lie between the equilateral and local shape. We also note 
that the sizes of the trispectra are even larger than those of the bispectra squared in this 
model. 

There are a lot of issues remaining to be investigated in the quasi-single field infiation 
models. For example, 

• Data analysis. Experimental constraints on non-Gaussianities depend on their detailed 
shapes and running. A variety of experimental methods have been applied to the local and 
equilateral bispectra [26143T] . It will be very interesting to constrain the family of new shapes 
that we find here, or to fit the parameter using the observational data. 

• Running of density perturbations. In this paper we have only considered the constant 
turn case. More realistically, we expect parameters to vary along the trajectory. As we noted, 
such a non-constant turn results in a running of the spectral index and non-Gaussianities. 
It is worth to investigate this running effect in more details. 
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• Trispectra. In this paper, we obtained the order of magnitude estimate of the trispectra 
in quasi-single field inflation. However, the zz-dependent coefficients in the trispectra, and 
more importantly the shape of the trispectra still remain to be calculated. 

• Multiple isocurvature directions. In this paper, we considered one massive isocurvaton. 
If there exist multiple massive isocurvature directions, instantaneously along the infiaton 
trajectory, one can find a two-dimensional hyper-surface where there is only one effective 
isocurvaton. If this hyper-surface is not changing with time, the model belongs to the 
two-field model such as the one we considered here. If the hyper-surface is changing with 
time, the situation becomes more complicated. It is worth to investigate the observational 
consequences of such a case. 

• Other quasi-single field mflaton models. The turning trajectory model that we studied 
is one simple example of the quasi-single field inflation models. The coupling between the 
infiaton and isocurvatons can be introduced through the kinetic terms in a more general way 
which may or may not be described by turning trajectories]^ Such couplings can even be 
introduced through other types of couplings that do not involve the kinetic terms. Which 
correlation functions/non-Gaussianities are enhanced is determined by the structure of these 
couplings. 

• String cosmology. As we discussed in the Introduction, quasi-single field inflation is a 
natural picture for inflation in string theory and supergravity. Fields with mass of order H 
are ubiquitous and in fact are a common hazard to inflation model building. We have seen 
that, as isocurvatons, such fields can have important consequences on density perturbations. 
It thus becomes very interesting to build explicit quasi-single field inflation models from 
string theory, in terms of either turning trajectories or more general couplings. Because the 
shapes of the non-Gaussianities depend very sensitively on the mass, and the sizes depend 
on the strength of the couplings, we have the opportunities to probe such fields through 
experiments. 

• Generalizations. As the number of fields increases, it becomes a logical possibility that 
we can have more than one infiaton fields. It is worth to generalize the formalism developed 
in the current work to these more general multiple field models. It is also worth to apply 
the perturbative method used here to multifield inflation models. 
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A The full Lagrangian up to third order 

In this section, we derive the full third order action in two different gauges mentioned in 
Sec. 12.21 We also show that (I2.13p is the leading order interaction. We set Mp = 1 in this 
Appendix. 

A.l Spatially flat gauge 

The full action is 



S Sg -\- Si 



(A.l) 



where 




(A.2) 



and 




(A.3) 



which is given in (12. ip . Using the ADM metric. 



ds"^ = -N^dt^ + hij{dx' + N'dt){dx^ + N^dt) , 



(A.4) 



the action becomes 




where the index of A^* can be lowered by the 3d metric hij and i?*^^-* is the 3d Ricci scalar 
constructed from /i,,. The definition of Ej~ and E are 



Eij = -{hij-ViN,-V,Ni) , 
E = Eijh'^ . (A.6) 
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We choose the following spatially flat gauge: 

hij = a\t)6ij , 9 = 9o{t) + 69 , a = a^it) + 5a . (A.7) 

For the constant turn case o"o(t) = const.. 

The constraint equations for the Lagrangian multiplers and Ni are 

f)C 1 

Ri^) +2C^ + 2N^-—{E,,E'= -E') = 0, (A.8) 

dC 

V^[N-\E'^ -h'^E)\+N^ = 0. (A.9) 

In the ADM formalism, to expand the action to the third order in perturbations, it is 
sufficient to solve the Lagrangian multipliers A^ and Ni to the first order in perturbations. 
So we expand 

N = l + ai, Ni = diiJi + N^^\ (A.IO) 

where diN-'^^ = 0, and plug them into (lA.SP and flA.9p . The solutions with proper boundary 
conditions are 

2 2 

= (gH^ - R'Ol) "1 ~ ^ (^R^9o69 + R9l5a + V^^59 + V'5a) . (A. 11) 

We then plug these solutions into the action and expand up to the third order in perturba- 
tions. The first order terms give the equation of motion for 6'o(t) and cro(^)- The quadratic 
order terms are 

C, =yR'69' - "-R^dM? - - (3e - + ^V)H') R"^(^" 

+ 2a^R9oS95G - 2ea^R9oH596a . (A. 12) 
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The third order terms are 

+ + 4m^^^ + ^^^^^ 

R IRQ 1 

-^5adi69di69 :T—5adi59diip 7rdiipdi5a5a 

o/' 



0,^ 

where 

A. 2 Uniform inflaton gauge 

In the uniform inflaton gauge, 

0{k, t) = 0o{t) , (7(x, t)^ao + 5(7(x, t) , (A.15) 

/ii,(x,t) = a2(t)e2«("'*)5i, . (A. 16) 

Again solving the constraint equations to the hnear order in perturbations, we get 

= A^f^ = 0, V'^-l + X, (A.17) 

where 

d\ = a'6(C - ^Sa) . (A.18) 

Plug these into the Lagrangian and expand. The linear order terms are consistent with the 
Hubble equation and give the equation of motion for cro(t). For the quadratic and cubic 
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order perturbations, after integrations by parts, we get 



H 



(A.19) 



where 



6 



c 



c 



-(^iH (t7 - C) + a6adi6a{^ - d,x) + Scr (-f- + 3C) 



H 



H 



H 



R R R 



2 4a 



5C 



4^ H Aa?H'^ 



1 



[d~-d.d,{d.CdjO-{d,Cf] 
d-\d,djCdidjx-d'Cd\) , 



2a( ^^^x + ifS^X-eS^C 



and X is defined in (lA.lSp . 



(A.20) 



(A.21) 
(A.22) 



A. 3 Estimate all cubic terms 

As we have mentioned in the Introduction, the contribution of the 6a^ term to f^L can be 
estimated as follows. V" / H and 9q/H are the dimensionless coupling constants for the cubic 
isocurvature interaction and the transfer vertex, respectively. Note that C ~ a/^- Then the 
three-point function is 

3 



V" I On 



H \H 

According to the definition f l4.15p . 



p; 



,3/2 



1 V" 



9n 



c 



1 V' 



(A.23) 



NL 



H \ H 



(A.24) 
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Other terms can be estimated similarly. For example, the 59' term in the spatially flat 
gauge contributes 0{e^)P^^'^ to (C^), where 0{e^) collectively denote second order in slow-roll 
parameters. So the 59' terms leads to fjsii ~ C^(e^)) which is negligible. In Table [1] and |2l we 
summarize the orders of magnitude of /tvl contributed from each term in the two gauges. 



term 
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59' 
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Table 1: The order of magnitude estimation for the contributions to Jnl from each term in 
the Lagrangian. The first row lists terms in the spatially fiat gauge. The second row lists 
the order of magnitude for /tvl contributed from the corresponding term. In this table, e 
collectively denotes all slow roll parameters. 
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Table 2: The same table for the uniform infiaton gauge. In the last entry, the 1st term is 
due to the field redefinition in 3pt, the 2nd is due to an extra cubic term from the field 
redefinition in the transfer vertex. 

As we can see the only term that can make /atl ^ 1 is the 5a' term. The rest of terms 
are corrections suppressed by either the slow- roll parameters or P^^"^ ~ ^l\f^ ~ H^/{R9o)- 

New cubic terms will arise after we convert the Lagrangian density to the interaction 
Hamiltonian density. Here we show that they do not change the order of magnitude estimate 
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performed above. The Lagrangian density in our case is of the form 



A = /o«' + /o/3' + /2 , 

£3 = hia^ + h2a + jj^ + j2$ + hs , (A.25) 

where a and P represent ( and Sa respectively, and the subscripts in various functions denote 
the orders of a and /3 in these functions. Following the same description in Sec. |2l we get 
the following interaction Hamiltonian density, 

^2 = /o«/ + /o/3/-/2 + |^ , 
= -giai , 

= -h,aj-(h2-^)ai-jj]-j2$i-h + ^ . (A.26) 
Jo ^/o 



Comparing to ( lA.25p . the new terms appeared in is 

foaj , (A.27) 



9! fSC2^' 



9ihi . 


(6C2 


Jo 


U2 


9ih2 


f6C2 


2/0 





2/0 V 

and in Ti^ are 

hiaj , (A.28) 
h2ai . (A.29) 
As long as 6C2/C2 ^ 1, these new terms will not be more important than the original ones. 

B Details of terms in various forms 

In this appendix we give the details of the terms in the factorized, commutator and mixed 
forms. 

B.l All 10 terms in the factorized form 

The three-point function is the sum of the following ten terms, 

/TT\3 10 

For simplicity, we only indicate the conformal time variables, such as ri or fi, once in each 
integrand. It applies to each factor before it. 
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rfri a^C2Vp^u'*{Ti) I dT2 a^C2Vp^u'*{T2) 



T3 



rfra a C3V*^v*2Up3(r3) / a C2vl^u'*^{u) 



P3 P3 ^ 



X (27r)^(5^(^pi) + 5perm. 



(B.9) 



(B.IO) 



(10) = -12up^Up^Up^{0) 

- /-O 

X Re 



/u /-n 
(iri a^C2VpiMp* (ri) / dT2 a? C2Vp^u'p^{T2) 
-00 J —00 

dT3 a^C2Vpy*^{T3) / dn a'^C3v;^v*yp^{r4) 

■00 J —00 

X (27r)^5^(^ pj) + 5 perm. 



(B.ll) 
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B.2 All 3 terms in the commutator form 

For convenience, we also list all three terms in the commutator form here. 



-12{ S-] Up^Up^Up,{0) 



Re 



/O pTl rT2 rT3 ^ 

dTi / dT2 / dTs / dT4 Y\ {a^{Ti)c2(Ti)) 
-OO J —OO J —OO J —OO 



/ ^C'iij2) . , nC3(t"3)j5 , xC3(t-4) 
0(^2) — 7 — + 0(^3) — - — -B + a(T4) ^ — rC- 



C2(,r2j ' ' C2(r3) 

X (27r)^(5^(^ pj) + 5 perm. , 



C2(r4) 



(B.12) 



where 



A = 



B = 



C 



Ki('^l) - C-C.) (^^pi(ti)^;;^(t2) - C.C.) (^^p3('^2)^^;3(t4X3(t4) - C.C.) 

Vp2{r2)vl^{n)u'*^{T3) , 

Ki('ri) - C.C.) {Up^{T2) - C.C.) {v*^{ri)v*^{T2)Vp,{T3)Vp^{T3) - C.C.) 
t^Ps{n)Vp^{T4)Up^{u) , 

- Ki('ri) - C.C.) {Up^{T2) - C.C.) {u'p^{t3) - C.C.) 

^Pi ('^iX ("^sX (rs)^;^! (t4)t;p2 (t4)t;p3 (T4) . 



(B.13) 
(B.14) 
(B.15) 



B.3 All terms in the mixed form 

We introduce a cutoff Tc to separate the IR and UV region. A quartic integral becomes 

/O /"Ti /•T2 /-ra 

dn / dT2 / dT3 / dT4 f{Ti,T2,T3,T4,) 
-00 J —00 J —00 J —00 

rO PTl pT2 rT3 

= dn dT2 / rfT3 / rfT4 /(Ti,r2,T3,T4) 

i/Tc «^Tc «^Tc »/Tc 

/(ti,T2,T3,T4) 



+ 
+ 





rT2 


/ c?r2 i 




/tc J 




l"ri 




/ rfr2 


/ dT3 




J —00 



dTi 

dTi 



J Tc -J Tc L-J —00 •/—CO 

rO r /-Tc PT2 i-Tz 

+ / dTi I dT2 I dT3 I dT 

J Tc L^— 00 J —00 J —00 



f{ri,T2,T3,T4) 
/(ti,T2,T3,T4) 



+ 



n pT2 pT3 

dn I dT2 / dT3 / ciT4 /(ti,T2,T3,T4) . 

00 J —00 J —00 J —00 



(B.16) 



We label the above five terms as (a), (b), (c), (d) and (e), respectively. 

We now write the three-point function in a mixed form. In the IR region, we write them 
in terms of the commutator form, so that the cancellation of the leading terms, especially the 
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cancellation of the IR spurious divergences, is manifest; in the UV region, we write them in 
terms of the factorized form, so that in each factor the fast UV convergence can be achieved 
by a Wick rotation and there is no spurious divergence for any momentum configuration. 

We start with the three terms in the commutator form flB.13l) - flB.15l) . We expand the 
expressions at a particular layer when the factorized form is need. To factorize the integral, 
combinations of different terms from all three terms. A, B and C, including the complex 
conjugates and the momentum permutations, are needed. It is always possible to make the 
required part of the integral factorize to a level indicated by the ten terms in Appendix IB.lt 
which cannot be further factorized. Hence the UV behavior becomes as good as those 10 
terms. At the same time, the IR convergence is as good as the three terms in Appendix IB. 2[ 

For the (a)-term in (IB. 161) . we use the integrand in the commutator form. So we have 3 
terms, which we label as (al), (a2) and (a3). 

For the (6)-term, the separation is simple, we get two terms from (IB.13P and one each 
from ( iRTill and ( IBTTbI) . We label as (61), (62), (63) and (64). 

For the ic] term, we have 



(cl) = 12MpjMp2Mp3(0) 



X Re 



X 



(c2) 



X 



X 



X 



X 



dTi 

c 

(27r)353(^p,) + 5 perm 

i 



dT2 a^Ciin) a^C3(r2) (^^^(ri) - c.c.) (t^pi(T-i)w*^(r2) - c.c.) Vp^Vp.,{T2) 



dr^ a^C2V* u*jTs) a'^C2V* u* (t^) 



(B.17) 



dri 



Tl 



dT2 a^C2{Ti) a^C3{T2) {Up^in) - c.c.) {Vp,{Ti)v*^{t2) - c.c.) Vp^V*^{T2) 

Tc -J Tc 

/Tc fTc 
dTs a^C2V*p^u'p^{Ti) / dr^ o? C2Vp^Up^{Ti) 
-oo J ~oo 

{2-Kf5^{Y^ p,) + 5 perm. , (B.18) 



which comes from (IB. 131) and we have used the complex conjugation to factorize (c2); 
(c3,4, 5) = l2up^Up^Up^{Q) 

pO PTl 

X Re I dri dT2 a^C2 {u'p^ - c.c.) Vp^in) a^C2 {u'p^ - c.c.) Vp^{T2) 



X 



oo 

Tc 



dT3 a C3Vp,Vp2Vp^,{T3) / dT4 a C2V* u'*{t4) 



dn 

dT3 



T3 



dTi a^C3Vp^Vp^v*p.^{n) a^C2Vp^u' {n) 



oo 
T3 



dTA a^C2V*u' (ts) a^C3Vp^Vp^Vp.^{T4^) 
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X (27r)353(^pi) + 5 perm. , (B.19) 

i 

where the contributions from ( IB. 141) and ( IB. 151) are combined to factorize the 1st term above, 
and the 2nd and 3rd terms (which cannot be further factorized) come from ( IB. 141) and (IB. 151) . 
respectively. 

For the ((i)-term, we have 

((il — 6) = 12up^Up^Up.j^{Q) 



X Re 





dri o?C2 {u'p^ - c.c.) Vp^{Ti) 



Tc 

Tc PTc rT3 



dT2 a^C2Vp^u' {t2) / dr^ I dr^ a^c^v* v* Vp^{T^) a^C2V* u'* (u) 



X 



OO J —OO J —OO 

Tc PTc PTg 

dT2 a^C2Vp.^u'p^{T2) / dr^ I dr^ C2Vp^u*^{T3) a'^OiV*_^v*^v*.^{n) 

■OO J —OO J —OO 

/Tc PTc PTS 

dT2 a'^C3V*^V*^V*^{T2) / rfrg / dT4 tt^ C2Vp2u'p^{T3) C2Vp.,u'p^{Ti) 
-OO J —OO J —CO 

dT2 dT3 dr^ a^C3V*^Vp^Vp^{T2) a^C2V*y*^{T3) a^C2V*^u'*^{T4) 

-OO J —OO J —OO 

Tc rT2 /"T3 

dT2 / dT3 / dn a'^C3V*,^Up^{T2) a'^C2Vp^Vp2V*^{T3) a^C2Vp^u'p.^{Ti) 

■OO J —OO J —OO 

/Tc l'T2 /"TS 

dT2 / dT3 / dr^ a^C2Vp2U*^{T2) C2Vp^u'*^{T3) a'^C3V*^V*^V*^{Ti) 
■OO J —OO J —OO 

{2TTf 6^{Y1 Pi) + 5 perm. , (B.20) 



where the contributions from (IB.lSp . (IB. 141) and (IB.15P are combined to factorize the first 
three terms above, and the last three terms (which cannot be further factorized) come from 
(iRTSj) . (IB141) and (IRTSD . respectively. 

For the (e)-term, we use the integrand in the factorized from in Appendix IB.li So we 
have 10 terms, which we label as (el), . . . , (elO). 



C Numerical integration and Wick rotation 

In the r — )■ — OO end of various integrations that we encounter in the in-in formalism, we 
slightly rotate r to the imaginary plane, r — )■ — oo(l ±ie), to achieve the convergence. If we 
have to compute the integration numerically, this is not easy to implement [32]. There are 
several ways to improve this situation. 

One way is to use the integration by part to increase the convergence speed of the 
integrand at the r — )■ — oo end Additional terms appeared in the integration by part 
are the integrand and its derivatives evaluated at a boundary. One may also choose a cutoff 
at some point tq and replace the integrand by an analytical approximation for t < tq. 
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In this paper, we use another way that is more efficient in terms of numerical integration. 
For an integral 

dn dT2--- dTnf{Ti,T2,---,Tn), (C.l) 



— oo 



we analytically continue — )■ ixi, so that the integral is equal t 

l>Xl rXn-l 

dxi I dx2--- / dxnf{ixi,ix2,--- ^ixn) ■ (C.2) 
±00 J ±00 J ±00 

After this Wick rotation, the original oscillating behavior of the integrands at the lower end 
becomes the exponential decay. Note that the lower limit for Xi is chosen to be either +00 
or —00 to make sure that the corresponding integrand is decaying as Xi — )■ — 00, and gives 
zero when evaluating the indefinite integration at Xi = —00. 

The reason that the above two expressions are equivalent is the following. We first look at 
the inner-most integral in (1C.1|) . Denoting the indefinite integration as F„_i(ri, ■ ■ ■ , r„_i, r„), 
we get F„_i(ri, ■ ■ ■ , r„_i, Tn-i) — F„_i(ri, ■ ■ ■ , r„_i, —00). For the type of integration that 
we have in the in-in formalism, as we take r„ — )■ —00, the lower end of -F„_i is quickly oscil- 
lating around a constant which we shift to zero in the definition of -F„_i. The slight rotation 
Tn — )■ — oo(l±ie) is chosen to make this oscillating piece go away, and the integration becomes 
F„_i(ri, ■ ■ ■ ,T„_i,T„_i). After the Wick rotation, such an oscillating behavior becomes an 
exponential behavior and is zero when we evaluate it at either 00 or —00. For the method to 
work, it is important that none of the other terms that originally goes to zero at infinity be- 
comes non-zero at infinity after Wick rotation. The upper limit F„_i(ixi, ■ • • , ix„_i) 
is the analytical continuation of -F„_i(ri, ■ ■ • , r„_i, r„_i). Performing the rest of the integra- 
tion in the same fashion, at the outer-most integral, because the upper limit is 0, we obtain 
Fi(0, ■ ■ ■ , 0) which is the same before and after the analytical continuation. 

For example, for ai > 02 > 0, 

fO 



dne-'"'^' / rfr2e*"^^2 = — . (C.3) 

02(01-02) 



After the analytical continuation, we have 



I / ' rfxse-"^^^ = —-^ . (C.4) 

J+00 02(0.1 — 0-2) 



00 



The integration (lC.4p is much easier to do numerically than (IC.Sp . 

To apply the above method to quasi-single field inflation and more general cases, some 
comments and discussions are in order here. 



^There are no singularities in the complex plane in our case. But it is an interesting question how this 
procedure will be modified if there are. 
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• The Wick rotation is not only useful to deal with the UV behavior of the integration, but 
also IR in some cases. If in the IR regime (r — )■ 0), there are both oscillation and exponential 
suppression in the integrand, caused by the mode function (12.261) . we can also perform a Wick 
rotation to find the total suppression factor. Instead of rotating 7r/2 in the complex plane, 
we can rotate another angle so that the integration after the Wick rotation only contains 
exponential suppression factor. This is useful when dealing with the m/H > 3/2 case. In 
this case, in the IR regime of the correlation functions, one needs to integrate 

adr (-r)-^"^+^ ■ ■ ■ oc [ dt e^*'^'^"^)^* ■ ■ ■ , (C.5) 



where a and b are some real numbers and we have approximated r oc e We rotate 

|zai/ — b\ 

Then the exponential factor in the integration becomes a suppression factor ~ ^-amt 
m ^ H. We also note that, when the modes are sub-horizon, they start to oscillate and 
their contributions are cancelled out. This corresponds to the scale t ~ 1/H. Combining 
both effects, we get factors of Boltzmann suppression e"™/^ for the correlation functions 
involving very massive modes. 

• If the upper limit of the outmost integral is non-zero, say Tc, we make a shift Tj — )■ Tj + Tc 
before the Wick rotation. This is used in the calculation of the mixed form. 

• This method is particularly easy if we know the analytical form of the integrand. This 
is the case in our paper. However, even if we do not, one can do the same Wick rotation for 
the equations of motion, such as (I2.22p and (I2.23p . 

-TY 1 ^ = , (C.7) 

dx^ X dx 

-TY - ^ + JT^^k = , (C.8) 

ax^ X dx ii^a;^ 

and solve the mode functions with the boundary condition 

Ruk , Vk — )■ -=xe^^ , as X — )■ — oo . (C.9) 

\/2k 



D Perturbation series in the in-in formalism 

In this appendix, we list various forms of in-in formalism. We start from [5J 



m)) = (oi 



Texp(^i dt'Hiit')^ Qj{t) Texp(^-i dt'Hiit')^ 



|0) . (D.l) 
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One can expand Eq. (ID.ip into series in two forms, namely, the factorized form and the 

commutator form. The nth order contribution of the factorized form is 

/■* _ rh _ rtn/2-i _ rt r^i /"*n/2-i 
^«(_l)"/2 dii dh--- din/2 dti dt2--- dtn/2 

J to J to J to Jto Jto J to 



X {Hiiin/2) ■ ■ ■ Hi{h)Qlit)Hj{h) ■ ■ ■ Hiitn/2)) 

/"* _ rh _ rtn/2-l-ni _ ft rh rtn/2-l + m 

+2Re^r(-l)"^+"/M rfti / dt2--- dtn/2-m dti dt2--- dt 

m=l •'til "'*n J to J to Jto 



X {Hi{tn/2~J ■ ■ ■ Hj{ti)Qlit)Hj{ti) ■ ■ ■ Hi{tn/2+m)) (D.2) 

for even n, and 

("±^/^ /■* _ /•*(n-l)/2-m _ ft ft(„-3)/2 + m 

2Re J2 / dh--- o?f(„+i)/2-m / dh--- d^n-l)/2+m 

™_T J ^ tn J tn J tn 



m=l 



X {Hi{^n+l)/2-m) ■ ■ ■ Hi{ti)Qi{t)Hi{ti) ■ ■ ■ Hji^n-m+rn)) (D.3) 

for odd n. 

The nth order contribution of the commutator form takes the form [6] 

ft ftl ftn-l 

/ dt, dt2--- dtn{[Hj{tn),[Hj{tn-i),--- ,[Hi{ti),Qi{t)]---]]) . (D.4) 

Jto J to Jto 

In the factorized form, all the integrands are under a unique integration. 

Each form has its computational advantages and disadvantages, as we discussed in detail 
in Sec. HI 

We can also rewrite all the integrals in terms of anti-time-ordered integral instead of 
time-ordered. Namely, the factorized form can be written as 



dtn/2 



pt pt pt pt pt pt 

,"(_l)-/2 / / ... dtr,/2 / dt, dt2--- 

Jto J ti Jtn/2-1 Jto J t\ Jtn/2~1 

X {Hj{h) ■ ■ ■ Hj{in/2)Ql{t)Hj{tn/2) ■ ■ ■ Hj{t,)) 

'"/2 pt nt pt pt pt pt 

+2Re5^z"(-l)"/2+"^ / dt, / dt2--- / din/2-m / dt, dt2--- dtn/2+m 

m=l *0 *1 -J tn/2-l-m ''to J tl '^<n/2-l + m 

X {Hi{t,) ■ ■ ■ Hi{tn/2-m)QAt)Hi{tn/2+m) " " " Hi{t,)) (D.5) 

for even n, and 

(n+l)/-^ pt pt pt pt 

-2Re ^ ,"(_l)-+("-l)/2 / dt,--- d\n+l)/2-m dt, - - - dt^n-l)/2+m 

m=l *0 >^*(n-l)/2-m -Jto >^t{n-3)/2 + m 

X {Hi{t,) - - - Hj{^n+l)/2-rr.)Qlit)Hj{t(^n-,y2+n.) " " " Hi{t,)) (D.6) 
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for odd n. And the commutator form can be written as 




dtn{[Hi{h),[Hj{h), - ■ ■ ,[Hi{tn),Qi{t)] - ■ ■]]) . (D.7) 
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